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Abstract 

The optical conductivity of graphene, bilayer graphene, and graphite in quantizing magnetic 
fields is studied. Both dynamical conductivities, longitudinal and Hall's, are analytically evaluated. 
The conductivity peaks are explained in terms of electron transitions. We have shown that trigonal 
warping can be considered within the perturbation theory for strong magnetic fields larger than 
1 T and in the semiclassical approach for weak fields when the Fermi energy is much larger than 
the cyclotron frequency. The main optical transitions obey the selection rule with An = 1 for the 
Landau number n, however the An = 2 transitions due to the trigonal warping are also possible. 
The Faraday/Kerr rotation and light transmission/reflection in the quantizing magnetic fields are 
calculated. Parameters of the Slonczewski-Weiss-McClure model are used in the fit taking into 
account the previous dHvA measurements and correcting some of them for the case of strong 
magnetic fields. 

PACS numbers: 71.70.Di, 78.20.Ls, 78.67.Wj 
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I. INTRODUCTION 



Comprehensive literature on the graphene family can be described in terms of the Dirac 
gapless fermions. According to this picture, there are two bands at the K hexagon vertexes 
of the Brillouin zone without any gap between them, and the electron dispersion can be 
considered as linear in the wide wave-vector region. For the dispersion linearity, this region 
should be small compared with the size of the Brillouin zone, i.e. less than 10 -8 cm^ 1 , 
providing the small carrier concentration n <C 10 16 cm -1 . Pristine graphene at zero temper- 
ature has no carriers, and the Fermi level should divide the conduction and valence bands. 
However, undoped graphene cannot be really obtained, and so far purest graphene contains 
about n ~ 10 9 cm -2 of carriers. Then the following problem appears — how do Coulomb 
electron-electron interactions renormalize the linear dispersion and does graphene become 
an insulator with a gap? 

Semiconductors with the gap are needed for electronic applications. Investigations of 
the graphene bilayer and multilayer are very popular as the gap appears when the bias is 
applied. We see how physics made a circle for the half of century returning to graphite 
studies. Here Slonczewski, Weiss, and McClure (SWMC) should be mentioned because they 
have stated the description of a layered matter- with interactions strong in the layer and 
weak between layers. 

The most accurate investigation of the band structure of metals and semiconductors is 
a study of the Landau levels through experiments such as magneto-optics^— and magneto- 
transportii"— . In magnetic fields, the classical and quantum Hall effects are observed, as well 
as the polarization rotation for transmitted (the Faraday rotation) or reflected lights (the 
Kerr rotation). However, the interpretation of the experimental results involves a significant 
degree of uncertainty, because it is not clear how the resonances can be identified and which 
electron transitions they correspond to. 

The theoretical solution for the band problem in magnetic fields often cannot be exactly 
found. A typical example is presented by graphene layers. For bilayer graphene and graphite, 
the effective Hamiltonian is a 4 x 4 matrix giving four energy bands. Fig. Q] shows the 
nearest two bands of the level structure together with the semiclassical orbits. The trigonal 
warping described by the effective Hamiltonian with a relatively small parameter 73 provides 
an evident effect (see the right panel). Another important parameter is the gate-tunable 
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Figure 1: (Color online) (a) The energy dispersion e(k,a) of two nearest bands (the electron band 
shown in solid line and the hole band in dashed line) in bilayer graphene for two polar angles a 
with the local extrema at k ^ ("mexican hat") represented. The band parameters are given in 
the figure, others are 70 = 3.05 eV, 71 = 360 meV, 74 = —150 meV. (b) Cross-sections k(a,e) of 
the electron band for energies of 80 meV (dashed-dotted line) and 40 meV (solid line). 

bandgap U in bilayer graphene. In this situation, the quantization problem cannot be solved 
within a rigorous method. To overcome this difficulty several methods have been proposed 
for approximate^^—, numerical^ - —, and semiclassical quantization^^. 

The present paper is organized as follows. In Sec. II we recall the electron dispersion 
in the graphene, bilayer graphene, and graphite. In Sec. Ill the optical conductivity and 
light transmission are discussed. In Sees. IV and V we describe in detail the quantization 
in magnetic fields. In Sec. VI the longitudinal and Hall conductivities as well as the 
Faraday/Kerr rotation are described. 

II. ELECTRON DISPERSION IN GRAPHENE FAMILY 
A. Electron dispersion in graphene 

The symmetry of K point is C^ v with the threefold axis and reflection planes. This group 
has twofold representation with the basis functions transforming each in other under reflec- 
tions and obtaining the factors exp (±27rz/3) in rotations. The linear momentum variations 
from the K point p± = ~^ip x — p y transform in a similar way. The effective Hamiltonian is 
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invariant under the group transformations, and we have the unique possibility to construct 
the invariant Hamiltonian linear in the momentum as 

ff(P>=( ° vp+ ), a) 

where v is a constant of the velocity units. The same Hamiltonian was written using the 
tight-binding model. 

The eigenvalues of this matrix give two bands 



£i,2 = Tv^jvl+Vl = Tvp, 

where the sign =p corresponds to holes and electrons. The gapless linear spectrum arises as 
a consequence of the symmetry, and the chemical potential at zero temperatures coincides 
with the band crossing due to the carbon valence. The cyclotron mass has the form 

1 dS(e) e 

m £ = o j — = ~2 ' 

2n de v 2 

and the carrier concentration at zero temperature n(fi) = fj 2 /nh 2 v 2 is expressed in terms of 
the chemical potential //. 

Tuning the gate voltage, the linearity of the spectrum has been examined in the 
Schubnikov-de Haas studies 2 ^ with the help of the connection between the effective mass 



and the carrier concentration at the Fermi level m(/i)t> = ^hy/im(n). The "constant" pa- 
rameter v was found to be no longer constant, but at low carrier concentrations n ~ 10 9 
cm -2 , it exceeds its usual value v = 1.05 ±0.1 x 10 8 cm/s (at concentrations n > 10 11 cm -2 ) 
by the factor of 3. 

This is a result of electron-electron interactions which become stronger at low carrier 
concentrations. The logarithmic renormalization of the velocity was found by Abrikosov 
and Beneslavsky in Ref.— for the 3d case and in Refs.— >^ for 2d graphene. Notice, that no 
phase transition was revealed even at lowest carrier concentration. We can conclude that 
the Coulomb interactions do not create any gap in the spectrum. 



B. Electron dispersion in bilayer graphene and graphite 

Bilayer graphene has attracted much interest partly due to the opening of a tunable gap in 
its electronic spectrum by an external electrostatic field. Such a phenomenon was predicted 
in Refs.— ^ and was observed in optical studies controlled by applying a gate voltage^— . 



4 



The graphene bilayer lattice is shown in Fig. [2j Atoms in one layer, i. e., a and b in the 
unit cell, are connected by solid lines, and in the other layer, e. g., ai and bi, by the dashed 
lines. The atom a (a x ) differs from b (b x ) because it has a neighbor in the adjacent layer, 
whereas the atom b (bi) does not. 




H(p) 



(2) 



Figure 2: Bilayer lattice 

The effective Hamiltonian of the SWMC theory can be written^ 2 .^ near the KH line in 
graphite as 

/ 75 vp+ 7i 74^P-/7o\ 

vp- 72 lAvp-ho 73^P+/7o 

7i 74^P+/7o 75 VP- 
\liVp+/^ 7 3 t;p_/7o vp+ 72 J 

where p± = ~^ip x — p y are the momentum components and jj are the functions of the p z 
momentum in the major axis direction, 

72 = 272 cos (2p z d ) , 75 = 275 cos (2p z d ) + A , 
7i = 27j cos (p z d ) for 2 = 1,3,4, 

with the distance do = 3.35 A between layers in graphite. The nearest-neighbor hopping 
integral 70 ~ 3 eV corresponds with the velocity parameter v = 1.5ao7o = 10 6 m/s and 
the in-layer inter atomic distance a = 1.415 A. The Hamiltonian (|2j) is represented in a 
somewhat different form than that in Ref.-. The relations between the hopping integrals 
in these forms are given in Table [B The recent estimate^ 1 ^ of the parameters agrees with 
those given in the Table. 

The electron spectrum of graphite is shown in Fig. [3J There are four levels labeled by the 
number s = 1,2,3,4 from below at any momentum. As a consequence of axial symmetry, 
twofold degeneration e 2 = £3 exists at p x = p y = 0, i. e., on the KH line. 



Table I: The parameters of the Hamiltonian, Eq. 
obtained in the experimental works, all in meV. 



their values in the SWMC model, and 
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a SWMC,Ref. 1 , "Mendez et al, Ref£, c Doezema et al, RefA d Charlier et al, Ref^i 
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Figure 3: Band structure of graphite 



In bilayer graphene, every layer has only one neighboring layer. Therefore, we have to 
set 72 = 75 = and to substitute 7, = 7, for i — 1, 3, 4 in Hamiltonian (T2J). The parameter 
U can be included also in the bilayer Hamiltonian as a result of the gate voltage. Then, the 
gap appears between the 82 and £3, and these bands acquire the form of "mexican hat" . An 
important point is that two points, K and K', are in the Brillouin zone transforming each in 
other under reflection. Such a reflection changes the U sign giving two different dispersion 
laws at the K and K' points. 
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III. OPTICAL CONDUCTIVITY 

We use the general expression for the conductivity as a function of the electric field 
frequency u and wave vector k in the form^^ 



Mp+) - e m (p_)]{(w + i5) 2 - [e n (p+) - £ m (p-)] 2 } J ' 

valid in the collisionless limit (u, kv) ^> r , where r -1 is the electron relaxation frequency, 
P-t = p ± k/2, and v l mn is the matrix element of the velocity operator 



determined by Hamiltonians (fT]) or ([2]). Hitherto, we did not use any peculiarities of the 
graphene spectrum The expression acquired only the factor 4 due to summation over spin 
and over six points of the K type (two per the Brillouin zone). 

The first term in Eq. ([3]) corresponds to the intraband electron-photon scattering pro- 
cesses. In the limit of the high carriers concentration (T, Ep) 3> kv , it coincides with the 
usual Drude-Boltzmann conductivity, if the substitution uj — > u + ir^ 1 is made. The second 
term owes its origin to the interband n — > m transitions with the infinitesimal 5 determining 
the bypass around the pole while integrating over the momentum p. The real part of this 
contribution is reduced to the well-known expression for the absorbed energy due to direct 
interband transitions. 

A. Optical conductivity of graphene 

For optical frequencies w > to, we can integrate in Eq. (JH]) over the angle and write the 
conductivity as 



(Tij(ui, k) 




(3) 



v = dH(p)/dp 



(4) 





— oo 







using the variable e = vp. 
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The intraband term can be integrated once more, 



° mtra {") = fir 2 'f r -n ln ( 2 cosh > ( 6 ) 
7m(w + ?r l ) 21 



where we write u + ir 1 instead of a; to take the small relaxation frequency into account. 
This Drude-Boltzmann conductivity at low temperatures T <C /x takes the form 

<j intra (u) = n . (7) 

In the opposite limit of high temperatures, the intraband conductivity Eq. (J6]) becomes 

2?> 2 T In 2 

= 2 ; 2 . (8) 

The temperature dependence of the relaxation rate in graphene is discussed theoretically in 
RefA 

The interband contribution in Eq. (j5]) integrated at zero temperatures gives 



inter/ \ 6 \ n( n \ 1 1 + 2 /i) 

a r (u) = — 9{u)-2a) In 

Ah ,) - 



(9) 



2tt (w - 2/z) 2 

where the function expresses the threshold behavior of interband electron transitions at 
uj = 2fi. The temperature smooths out all the singularities 

9(u - 2fjt) -»> - + - arctan - ~ 2// (10) 
2 7r 2T 

(w - 2/i) 2 ->■ (w - 2/x) 2 + (2T) 2 . 

The main issue should be emphasized. In high frequency region u ^> (T, /i), the iterband 
transitions make the main contribution into conductivity 

a{u) = h 

having the universal character independent of any material parameters. This frequency 
region is limited above by the band width of around 3 eV. Making use the universal con- 
ductivity, one can calculate the light transmission through graphene^ in the approximation 
linear in conductivity 

An e 2 

T = l Re <j(u) cos 9 = 1 — it— cos 9 , (11) 

c he 

where 9 is the incidence angle. The intensity of reflected light is quadratic in the fine 
structure constant a = e 2 /hc. In excellent agreement with the theory, for the wide optical 
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range, several experimental groups^— observe the light transmission through graphene as 
well as bilayer graphene where the difference from unity is twice as larger. It is exception- 
ally intriguing that the light transmission involves the fine structure constant of quantum 
electrodynamics having really no relations to the graphene physics. 

For graphite, the value of <Jd = e 2 /4hdo plays the role of the universal dynamical 
conductivity, where do is the distance between layers. As shown experimentally^ and 
theoretically—, the dynamical conductivity of graphite is close to this universal value in 
the frequency range 0.1-1 eV having the kink singularity at the interband transition fre- 
quency lo = 271. 



IV. GRAPHENE IN MAGNETIC FIELDS 



In the presence of the magnetic field B, the momentum projections p + and p- become 
the operators with the commutation rule {p + ,p-} = —2ehB/c. We use the relations 



vp+ = ub a, vp- = ujb a 



involving the creation a + and annihilation a operators with ub = V\j2\e\hB/c . We will 
write only one of two x, y space coordinates including the corresponding degeneracy propor- 
tional to the magnetic field in the final results. 

For graphene, we search the eigenfunction of Hamiltonian ([1]) in the form 

«.<»> = ( C t?" ff • (12) 

where ip n (x) are orthonormal Hermitian functions with the Landau number n > 0. Canceling 
the Hermitian functions from the equations, we obtain a system of linear equations for the 
eigenvector C sn 

I x { ™ =0 (13) 
giving the eigenvalues 

e sn = TUB\/n 

with s — 1, 2 and n — 0, 1, 2... For n — 0, there is only one level £10 = with Cq = 0, Cq = 1 
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as follows from Eq. (jl~2l) . The wave function columns write 



C - 1 j 1 ^ 1 HA\ 

= —= < and (14) 

I -1 1 



for s = 1 and s = 2 and n — 1,2.... 



V. GRAPHENE LAYERS WITH TRIGONAL WARPING IN MAGNETIC 
FIELDS 



We search the eigenfunction of Hamiltonian ([2]) as a column 



CsnPn(x) 



(15) 



x 



One sees the every row in Hamiltonian fl2]) becomes proportional to the definite Hermitian 
function if the terms with 73 are omitted. We will show that the terms proportional to 73/70 
can be considered within the perturbation theory or the semiclassical approximation. 

Canceling the Hermitian functions from the equations, we obtain a system of the linear 
equations for the eigenvector C sn 



75 - £ 
7i 



72 -E 

\ u^\/n — 1 



7i 



n 

u^y/n 

75 - e u B y/n - 1 



C 



1 

ST* 



X < 



r<2 







(16) 



C 



i 

sn 



u B vn - 1 72 - £ ) 

where the band number s = 1,2,3,4 numerates the solutions at given n from the bottom, 
ug = vyj2\e\%B /c and U4 = ^ a uj b /^q. 

The eigenvalues of the matrix in Eq. f|T6|) are easily found, they are shown in Fig. H] as a 
function of the momentum p z . For each Landau number n > 2 and momentum p z , there are 
four eigenvalues e s (n) and four corresponding eigenvectors, Eq. (|15|) . marked by the band 
subscript s. We use the notation \sn) for levels. In addition, there are four levels. One of 
them is 

ei (n = 0) = 72 (17) 
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Figure 4: (Color online) Landau levels in graphite e sn for n from to 4 in four bands s =1,2,3, and 
4 (in dotted, solid, dashed, and dash-dotted lines, correspondingly) as functions of wave vector k z 
along the KH line in the Brillouin zone (K = 0, H = 7r/2do) in the magnetic field B = 7 T with 
the SWMC model parameters given in Table l.The main electron transitions shown in the right 
panel below 100 meV occur between the levels with the selection rule An = 1 , see text. 



for n = with the eigenvector Co = (0, 1, 0, 0) as is evident from Eq. (fl~5l) . It intersects the 
Fermi level and belongs to the electron (hole) band near the K (H) point. The other three 
levels indicated with n = 1 and s = 1, 2, 3 are determined by the first three equations of the 
system ([16]) with = 0. 

The |21) level is close to the 1 10) level. In the region of p z , 71/ cos2p 2 <io 3> 72, where the 
electrons are located, this level has the energy 

e 2 (n = l) =72-2^1. 

7i7o 

In the same region, the two closest bands (s = 2, 3) with n > 2 are written as 



£2» 

within accuracy of (74/70) 2 . 



72 



w|7 4 



UJ 



(2n-l)T—Vn(n - 1) 

7i7o 7i 



(18) 
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A. Perturbation theory for matrix Hamiltonian 



Due to a double degeneracy existing on the KH line, the effect of the trigonal warping 
becomes essential. A simplest way^ to evaluate the corrections resulting from the warping 
73 consists in the consideration of the Green's function having the poles at the electron 
levels. 

The Green's function of the unperturbed Hamiltonian writes using the functions in Eq. 
(USD as 

Gf(e,x,x') = ^ niX) _ rIiX,) - (19) 

£ £ sn 

sn 

The corrections to the levels can be found in the iterations 

G m+1 (x,x') = J d 2 x"G (x,x")V(x")G m (x",x') , (20) 

where V(x) has only two matrix elements V 42 = ujb13Cl + /lo and V 24 = V 42 * in the Hamil- 
tonian (j2J). 

In the second iteration, we get the corrections 

d 2 Xl d 2 x 2 Gf(x, x 1 )V 42 (x 1 )Gf(x 1 , x 2 )V 24 {x 2 )Gf{x 2) x') 

and the similar term with the superscript substitution 2 H 4. The matrix elements of the 
perturbation V are easily calculated with respect to the Hermitian functions in Eqs. f fT9|) 
and (EES}, and we obtain 



V (n - 2)\Ct n C^ n _^ n {xW/ n {x') 
7o / (e - e sn )(e - e s/}n ^s)(e - e sn ) 

for the diagram shown in the upper part of Fig. This correction plays an important role 
near the poles of the Green's function. For this reason, for e close to e sn , the e value in the 
second factor of the denominator can be replaced by e sn . Thus, the total Green's function 
(with the correction) has the structure 

1 | 5 

£ £ sn \~ ^sre) 

which can be rewritten to the second-order terms in 5 as 

1 



£ ^ sn O 
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a 



Figure 5: Diagrams for the second iteration of the perturbation theory; corrections to the Green's 
function (a), corrections to the vertex in conductivity (b). 

Therefore, we can represent the correction as a shift Se sn of the poles (e — e sn — fe sri ) _1 



where the first term should be omitted for n — 3 < 0. In fact, our illustration is nothing 
but a calculation of the electron self-energy and the naive expansion of the denominator can 
indeed be replaced by summarizing of the corresponding diagrams. 
The corrected 1 10) level writes 



The 1 21) — level is very close to the level with n = 0, Eq. (123]) . 

Comparing the corrections, Eq. f l22|) . with the main contribution Eq. ffl8|) . we find, first, 
that the perturbation theory is valid when an expansion parameter (7371/70^5) 2 becomes 
small, i.e., for strong magnetic fields B > IT. Second, the effect of 74 is linear, whereas 
of 73 is quadratic in these constants. Therefore, the 74 constant is more essential for the 
electron levels in magnetic fields. 

Comparison shows that Eqs. (122]) and (T23l for levels give the same results as the numerical 
method of truncating the infinite-rink matrix in Ref.— . 

Note that the derived expressions are applicable as well to bilayer graphene while one 
includes the field U and substitutes 72 = 75 = and 7$ = ji for i — 1, 3, 4. In the simplest 
approach, when only main parameters 71 and U are holded, the magnetic levels e sn are 



with 




(22) 




(23) 
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determined by the equation 

[(U _ 6sn f - u l n ][(u + £sn ) 2 - uj 2 B (n + 1)] + ^(U 2 - e 2 J = . 

B. Berry phase, semiclassical quantization and Landau levels 

Alternatively, the semiclassical quantization can be applied for relatively weak magnetic 
fields when the cyclotron frequency is small compared to the Fermi energy. Then, we can 
use the Bohr-Zommerfeld condition as 

(24) 

Here S(e) is the cross-section area of the electron orbit in the p x .p y space for the energy e 
and the constant momentum projection p z on the magnetic field, n c is an integer supposed 
to be large. The integer T is the number of the smooth turning points on the electron 
orbit. There are two smooth turning points for the Landau levels and only one for skipping 
electrons reflected by the hard edge. 

We use the semiclassical approach for the magnetic field normal to the layered system 
when the in-layer momentum components p x and p y are only quantized and the size of the 
Fermi surface is small compared with the Brillouin zone size. Notice, that the 5(e)-phase 
depends on the energy. If the spin is neglected, 5 = and T = 2 for the Landau levels, 
and 5 = 1/2 and T = 2 for monolayer graphene. In these two cases, the semiclassical result 
coincides with the rigorous quantization and it is closely connected with the topological 
Berry phased. This 5-phase was evaluated for bismuth in Ref.—, preceding Berry's work by 
almost two decades, and it was considered again for bismuth in Ref.—. For graphite, the 
semiclassical quantization was applied in Ref.—. However, in the general case, the evaluation 
of the 5— phase is still attracted a widespread interest^— . 

The problem under consideration is described by the Hamiltonian in Eqs. (JTJ) or (J2]) 
rewritten in the form 

(V-p + r-e)^ = 0, (25) 

where p and V are the two-dimansional vector and matrix, correspondingly, with the in- 
layer components x and y. The column \1/ is labeled by the band subscript which we omit 
together with the matrix subscripts on T and V, summation over them is implied in Eq. 



ehB 



S(e) = 2tt 



r 

n c + — + 5(e) 
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(|25|) . Matrices T and V are the first two terms (of zero and first orders) in a series expansion 
of the Hamiltonian in the power of quasi- momentum p x and p y . 

In the magnetic field, the momentum operator p depends on the vector-potential A by 
means of the Peierls substitution, 

p = —ihV — eA/c, 

providing the gauge invariance of the theory. The magnetic field can also enter explicitly 
describing the magnetic interaction with a spin of particles. However, for the graphene 
family, the magnetic interaction is weak and omitted here. 

It is convenient to choose the vector-potential in the Landau gauge A x = —By, A y = 
A z = in such a way that the Hamiltonian does not depend on the x coordinate. We search 
the function \1/ in the form 

\1> = $exp (is/%) , 

where the function s is assumed to be common for all component of the column 
The function $ is expanding in series oiU/i: 




Collecting the terms with the same powers of h in Eq. ( 1251) . we have 

(V • p + T - s)<p m = -VV Vm _! . (26) 

For m = 0, we get a homogeneous system of algebraic equations for the wave function 
column if , 

(V-p + T-e) Vo = 0, (27) 
which has a solution under the condition 

Det(V-p + r-e) = 0. (28) 

This equation determines the classical electron orbit, e(p x ,p y ) = e, at the given electron 
energy e in presence of the magnetic field. On the other hand, the equation coincides 
with the dispersion equation since it does not contain the magnetic field. In 3d case, as in 
graphite, the dispersion depends also on the momentum projection p z on the magnetic field. 
Thus, our scheme does not requires the expansion in a power of p z . 
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The equations (126]) with m = 0, 1 give the wave function in the semiclassical 
approximation^. The quantization condition can be written as usual from the require- 
ment that the wave function has to be single-valued. Making the bypass in the complex 
plane around the turning points to obtain the decreasing solutions in the classically unac- 
cessible region, we obtain, first, T = 2 and, second, 5-phase as a contour integral along the 
classical orbit 

,( e ) = l Im /_^ W ^£, (29) 
2tt J (f* <f v y y dp x 



where v y = de(p x ,p y )/dp y . Using the Hamiltonian hermiticity, after the simple algebra (see 
Ref.— ), Eq. (|29l) can be rewrite in the gauge-invariant form 

dp 



5(e) = -*-Im / — - 



:<Po 



Vx — 
dp 



(30) 



where v = \Jv\ + v 2 and the integrand is called the Berry connection (or curvature). Ev- 
erywhere, the summation over the band subscript is implied. 

Now let us calculate the 5-phase for bilayer graphene. In the simplest case, omitting 73 
and 74, the effective Hamiltonian can be written as 

/ U q+ 71 \ 
q- U 
71 -U q- 
\ q + -U J 

where the parameter U describes the tunable gap as a result of the gate voltage and 71 is 
the interlayer nearest-neighbor hopping integral energy. The constant velocity parameter v 
is incorporated in the notation q± = vp±. The band structure is shown in Fig. The 



(31) 



minimal value of the upper energy 64 is \/U 2 + 7^, the £3 band takes the maximal value \U\ 
at q — 0. Here, the orbit is the circle defined by Eq. ff28|) . written in the form 



[(U + ef - ?\ [{U - ef - <?] - yi(e» -U') = 0. 
The eigenfunction ip of the Hamiltonian (131]) can be taken as 

/ (U-e)[(e + U) 2 -q 2 }\ 



(32) 



^0 



q 4 q 2 -(e + U) 2 ] 
li{U 2 -e 2 ) 
iiq+(u -e) 



(33) 
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with the norm squared 

<p* <p = l(e + Uf - q 2 ] 2 [(e - Uf + q 2 } 

+1 l(e-U) 2 [{e + U) 2 + q 2 ]. (34) 

The derivatives for Eq. (129]) are calculated along the trajectory where the energy e and, 
consequently, the trajectory radius q are constant. 




Figure 6: Band structure of bilayer graphene 

If the conditions \U\ < \e\ < \JU 2 + 7i are fulfilled, Eq. ( 132]) has only one solution for 
the radius squared 



q 2 = U 2 + e 2 + ^AU 2 e 2 + (e 2 - U 2 ^ 2 . 

The matrix V y = OH/ dp y in Eq. (129]) has four nonzero elements, V y 12 = V 21 = V^ 34 = V^ 43 = 
-1. 

Using Eqs. (152) and (15311. we find 

Im <p* V v ?p- = AUe(U - e) [{e + U) 2 - q 2 } . (35) 

This expression is constant on the trajectory as well as <Pq<Po, Eq. (1341) . Therefore, in order 
to find S, Eq. (129]) . we have to integrate along the trajectory 

rdp* 

J Vy 

This integral equals —dS(e)/de, where S(e) = nq 2 is the cross-section area, Eq. ( 124"]) . with 

dS(e) 2{q 2 + U 2 -e 2 ) +1 2 

~dT = n£ q 2 -U 2 -e 2 ■ (36) 
17 



Now we have to substitute Eqs. (I34I) - (|36|) into Eq. (129|) . Thus, we find the Berry phase 

— ell —ell 

6 W = g2 _ £ 2 _ ip = ^pp = p _ 

shown in Fig. [71 where <5-phase of bilayer graphene with trigonal warping is also shown; the 
detailed calculations will be published elsewhere. For the ungaped bilayer, U — 0, the Berry 
phase 5(e) = 0. The Berry phase depends on the energy and 5 = =Fl/2 at e = ±U. At the 
larger energy, e^> U, the Berry phase 5 — > =Fi7/7i- 

Substituting Eq. (1371) in the semiclassical quantization condition, Eq. ( IMl) . and solving 
the equation obtained for e, we get the energy levels as functions of the magnetic field. We 
have to notice that the Landau numbers n listed in Fig. [8] do not coincide with the numbers 
n c in the semiclassical condition ( 12^|) . The rigorous quantization shows that there are only 
one Landau level with n = and three Landau levels with n = 1—. These levels are not 
correctly described within the semiclassical approach. However, for n > 2, there are levels 
in all four bands s (two nearest bands with s = 2, 3 are shown in Fig. E]). They correspond 
with the semiclassical number n c = n — 1, and the semiclassical levels for the larger n are 
in excellent agreement with the levels obtained in the perturbation approximation. 



-0.05 




e (meV) 



Figure 7: (Color online) Semiclassical phase vs energy in the conduction band of bilayer graphene 
without trigonal warping (solid line) and with warping (dashed line). 
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VI. MAGNETO-OPTICS EFFECTS IN GRAPHENE LAYERS 



An important peculiarity of conductivities in presens of magnetic fields is an appearance 
of the Hall component a xy (u). The Hall conductivity violates the rotation symmetry of 
graphene around the major axis. This implies the rotation of the linear polarized electro- 
magnetic wave, i. e., the Faraday and Kerr effects for transmitted and reflected waves, 
correspondingly First of all, the electron transitions are possible between the levels with 
the neighboring Landau numbers n and various bands s, and therefore the resonance de- 
nominators A ss i n = e sn — e a ' n +i arise in the conductivity tensor. 

Calculations^ give the conductivities for graphite in the collisionless limit when the 
electron collision frequency T is much less than the level splitting 

\ tt/2 
<TxxM I . 4cj 2 BSr f Afss'n u ,2 

> =%o d -—y \ dz- — — \d ss , n \ 
i^M J * 2 ^ S ,J ^s'n (38) 

x [(u + iT + As^n)- 1 ±(u + iT- As^n)- 1 ] , 

where the integration is taken over the reduced Brillouin zone, < z < tt/2. Such integration 
is absent for graphene and bilayer. Here Af ss i n = f(e s i n+ i) — f(e sn ) is the difference of the 




Figure 8: (Color online) Energy levels e sn for the K valley in magnetic fields for bilayer graphene 
within the perturbation approximation (solid lines) and in the semiclassical approach (dashed- 
dotted lines); in the notation \sn), n is the Landau number and s = 1,2,3,4 is the band number, 
only two nearest bands (s = 2, 3) are shown at given n from to 7. There is only one level, |10), 
with n = and three levels (s = 1, 2, 3) with n = 1. The levels for the K' valley can be obtained 
by mirror reflection with respect to the e = axis. 
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Fermi functions and 

+(%/lo)(C sn C s , n+1 + Cg n Cg, n+1 ) 
is the dipole matrix element expressed in terms of wave functions (ITST) . These transitions 
are most intensive. They obey the the selection rule 

An = 1 , 

and will be referenced as the strong lines. The conductivity units here 

„2 



0d 



Ahd Q 

have a simple meaning, being the graphene universal conductivity e 2 /Ah multiplied by the 
number l/ofo of layers within the distance unit in the major axis direction. 

Besides, we have to take the renormalization of the dipole moments due to trigonal 
warping into account. This additional electron-photon vertex results in weak lines with the 
selection rule 

An = 2. 

We get this contribution by substituting 

instead of the matrix element in Eq. fl38|) and replacing the subscript n + 1 — > n + 2. We 
have to notice, that the 74 corrections give the linear (in small parameter 74/70) contribution 
to the conductivities at the main electron transitions with An = 1. The 73 corrections are 
quadratic, however, they result in an appearance of new resonant transitions with An = 2. 

There are also small so-called vertex corrections to the self-energy shown at the bottom 
of Fig. |5j They result from the quartet of the coupled Landau levels, which interfere while 
the selection rules An = 1 and An = 2 are allowed. 



A. Gapped bilayer graphene 

Graphene and bilayer graphene effect the transmission and the Faraday rotation in a 
linear order in the fine structure constant whereas the reflected light intensity is quadratic 
in a. Therefore, let us discuss the characteristics of the transmitted light through bilayer 
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Figure 9: (Color online) Transmission spectra of gapped bilayer graphene without and with trigonal 
warping (dashed and solid lines, correspondingly) at 10 T and U = 30 meV; the band parameters 
used are v = 1 x 10 8 cm/s, 71 = 360 meV, 74 = —150 meV, ef = 30 meV, others are listed in Fig. 
The relaxation frequency is supposed as V = 5 meV. 

graphene where the effects have a maximum value. For this case, Eq. (158]) is valid without 
the integration over the z momentum component. The conductivity units should be taken 
now as o"o = e 2 /4h. In the approximation linear in conductivities, the transmission coefficient 
T and the Faraday angle for the free standing bilayer write as 



Results of calculations are shown in Figs. |9] and [TOl The peaks in absorption, Fig. [9] 
correspond to the electron transitions. There is the series of seven lines in the 0.1-0.4 eV 
interval. They are doublets excited by the electron transitions of the type \2n) — » |3, n+ 1) 
and \3n) — > |2, n + 1) for n from 2 to 8. Two weaker lines at 350 and 380 meV are resulted 
from the 1 10} — > |31) and |21) — > 1 42) transitions, correspondingly. There is strongest line 
at 24 meV excited by the |21) — > 1 32) transition. All these lines obey the selection rule 



The very weak lines at 51 and 78 meV owe their appearance to the An = 2 transitions 
|21) |33) and |10) ->• |22). 

In general, the effect of the small constants 73 and 74 is seen more on the low levels 1 10) 
and |21). 




xy ■ 



(39) 



An 



1. 
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Figure 10: (Color online) Faraday rotation in gapped bilayer graphene; parameters used are the 
same as in Fig. [9j 
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Figure 11: (Color online) Real (a) and imaginary (b) parts of the longitudinal (xx, solid line) 
and Hall (xy, dashed line) dynamical conductivities calculated for one graphite layer in units of 
an = e e /4:h; Kerr angle (c) and reflectivity (d). The magnetic field B = 7 T, the temperature T = 
0.1 meV is less than the level broadening T = 3.5 meV. 



The transition frequencies in the Faraday rotation, Fig. [9]are determined by the derivative 
maximum values. 
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Figure 12: (Color online) Kerr angle and reflectivity at 10, 15, and 25 T. 
B. Graphite 

Using the conductivities Eqs. ( 138]) . one finds the complex bulk dielectric function = 
5ij + 47ii<Jij/u and the reflection coefficient and the Kerr rotation (see, e.g.,—) 

1 1 

R = 2^\ r +\ 2 + \ r -\ 2 ),®K = -arg(r_r|) , 

where r± — (1 — y/£±)/(l + y/e±) are the reflection Fresnel coefficients for two circular 
polarizations with e± = e xx ± e xy . 

The parameters of Eq. (j2J) used in the calculations are listed in Table fl] (see also Ref.— ). 
The hopping integrals 70 to 73 are close to the values determined in observations of the 
semiclassical ShdH effect. The Fermi energy equal to £p = —4.1 meV agrees at the zero 
magnetic field with the measurements of the extremal Fermi-surface cross sections and the 
masses of holes and electrons. Connections with the notation for similar parameters of the 
SWMC model are given in the "SWMC" line. The values of parameters 74, 75, and A 
determined in various experiments are very different; we use 75 and A obtained by Doezema 
et al4 (given in Tabled] in the "SWMC" notations) and take the approaching value for 74. 
In the quantum limit, when electrons and holes occupy only 1 10) and |21) levels, the Fermi 
energy must cross these close levels at the middle of the KH line. It means that the Fermi 
level becomes higher at such the magnetic fields taking the value Ef ~ — 1 meV. 

The results of calculations are represented in Figs. [TT¥T2l Let us emphasize that the 
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imaginary part of the dynamical conductivity is of the order of the real part. 

One can see in Fig. [TTl(a). that the averaged longitudinal conductivity calculated per one 
graphite layer tends to the graphene universal conductance. The main contribution in the 
sharp 16-meV line is resulted from the electron |21) — > 1 32) transition (15 meV) about the 
K point (see Fig. H]) where the 1 32) level coincides with the Fermi level (within an accuracy 
of the width V or temperature T). Then, the transitions 1 22) — > |21) produce the broad 
band. The low- frequency side of the band (23 meV, at the intersection of the |21) level with 
the Fermi level) contributes into the 16-meV line. In the same 16-meV line, the transitions 
1 32) — > 1 33) can contribute as well if the band 1 32) contains the electrons. 

The next doublet at 43 meV arises from the transitions 1 23) — > 1 32) and 1 22) — > 1 33) at 
the K point. The 68- meV doublet appears as the splitting of the 1 24) — > 1 33) (65 meV) 
and 1 23) — > 1 34) (69 meV) transitions due to the electron- hole asymmetry at the K point of 
the Brillouin zone. 

The 89-meV line is more complicated. First, there are the electron transitions 1 24) — > 1 35) 
(89 veV) and 1 25) ->■ |34) (90 meV) near the K point. Besides, the transitions |11) ->■ |10) 
(95 meV) near the H point make a contribution as well. All these lines obeying the selection 
rule An = 1 are strong. There are two weak lines in the frequency range. One ( 1 24) — > 1 32) ) 
is seen at 55 meV as a shoulder on the theoretical curve. Another, at 31 meV, results from 
the transitions 1 10) — > 1 32) near the K point. 

The positions of the lines for fields in the range of 10-30 T agree with observations of 
Refs.^. 

The optical Hall conductivity a xy (u) in the ac regime is shown in Figs. [H](a) andfTTl(b). 
The conductivities cr xx (u) and a xy (u>) allow calculating the Kerr rotation and the reflectivity 
as functions of frequency [see Figs. [11] (c) andfTTI (d)]. It is evident that the interpretation 
of the Kerr rotation governed by the conductivity a xy (u) is much more complicated in 
comparison with the longitudinal conductivity. The Kerr angle and reflectivity shown in 
Fig. [12] for the different magnetic fields demonstrate the strong field dependence of the 
magneto-optic phenomena. 
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VII. SUMMARY AND CONCLUSIONS 



In conclusions, we have evaluated the perturbation theory for the matrix Hamiltonian, 
which permits to calculate the corrections to eigenvalues resulting from the small matrix 
elements particularly from the trigonal warping. The trigonal warping in graphite can 
be considered within the perturbation theory at strong magnetic fields larger than 1 T 
approximately. For weak magnetic fields, when the Fermi energy much larger than the 
cyclotron frequency, the semiclassical quantization with the Berry phase included can be 
applied. We have found that the main electron transitions obey the selection rule An = 1 
for the Landau number n, however the An = 2 transitions due to the trigonal warping with 
the small probability are also essential. In graphite, the electron transitions at the K, H 
points as well as at intersections of the Landau levels with the Fermi level make contributions 
into conductivity. The good agreement between the calculations and the measured Kerr 
rotation and reflectivity in graphite in the quantizing magnetic fields is achieved. The SWMC 
parameters are used in the fit taking their values from the previous dHvA measurements 
and increasing the Fermi energy value for the case of the strong magnetic fields. 
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